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Abstract 

A theory of the competition between fusion and quasi-fission in a heavy fusing 
system is proposed, which is based on a master equation and the two-center shell 
model. Fusion and quasi-fission arise from a diffusion process in an ensemble of 
nuclear shapes, each of which evolves towards the thermal equilibrium. The master 
equation describes the diffusion of the nuclear shapes due to quantum and thermal 
fluctuations. Other crucial physical effects like dissipation, ground-state shell effects, 
diabatic effects and rotational effects are also incorporated into the theory. The 
fusing system moves in a dynamical (time-dependent) collective potential energy 
surface which is initially diabatic and gradually becomes adiabatic. The microscopic 
ingredients of the theory are obtained with a realistic two-center shell model based 
on Woods-Saxon potentials. Numerical calculations for several reactions leading to 
^^^No are discussed. Among other important conclusions, the results indicate that 
(i) the diabatic effects play a very important role in the onset of fusion hindrance 
for heavy systems, and (ii) very asymmetric reactions induced by closed shell nuclei 
seem to be the best suited to synthesize the heaviest compound nuclei. 
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1 Introduction 



The fusion of massive nuclei is currently a very hot topic in nuclear physics 
motivated by the search for new superheavy elements (SHE) [1,2,3]. Although 

Email address: alexis.diaz-torres@anu.edu.au (Alexis Diaz- Torres) . 



Preprint submitted to Elsevier Science 



9 February 2008 



the periodic table of the nuclei has been successfully extended during the past 
few years by the experimental discovery of new SHE [1,2,3], the understanding 
of the formation mechanism of these nuclei is still a challenge for theory. Very 
recent experiments [4] at LBNL in Berkeley (USA) with ^^Ca + ^^^U have 
not confirmed the production of the ^^^112 SHE reported by the Dubna group 
(Russia) [3]. 

In contrast to the fusion of light and medium mass systems, where the fusion 
cross sections are determined by overcoming or tunnelling through the external 
Coulomb barrier of the nucleus-nucleus potential {the capture process which is 
rather well understood within the coupled channels framework, e.g., see [5] and 
references therein), the fusion cross section for heavy systems can be inhibited 
by many orders of magnitude (fusion hindrance [6,7,8,9,10,11,12,13,14,15]) due 
to the competition between fusion and quasi-fission (reseparation before com- 
pound nucleus (CN) formation) after the capture stage. The understanding 
of this competition remains an unsolved problem and is important for opti- 
mization of the SHE production rate as well as for understanding the new 
experimental results on fusion of heavy nuclei [11,12,13,14,15]. In addition to 
the capture process and the subsequent evolution of the combined system from 
the contact configuration into the CN {the stage of CN formation), the cool- 
ing process of the excited CN by the emission of light particles and gamma 
rays in competition with fission {the deexcitation process) clearly affects the 
formation of the evaporation residues (ER) observed in the experiments. The 
survival against fission is optimized by minimizing the CN excitation energy 
and angular momentum. This process is broadly described within a statisti- 
cal decay model of atomic nuclei [16]. It is based on the assumption that all 
nuclear degrees of freedom are equilibrated after the formation of the CN and 
before it begins to decay. For an overview of the three stages of a heavy ion 
fusion reaction leading to the formation of a heavy ER, we refer to Ref. [17]. 

Most of the theoretical groups involved in SHE research have a similar view- 
point regarding the description of the capture stage and the deexcitation pro- 
cess, but there is no consensus [18] on the modelling of the intermediate stage 
of CN formation. The physics of the current theories sometimes contradict 
each other. For the sake of simphcity, most of them assume that the fusing 
system is in thermal equilibrium. Depending on the main coordinate for fusion, 
two sorts of models can be distinguished: 

• In the first type [19,20,21,22,23,24,25,26,27,28], the fusion occurs in prac- 
tice essentially along the radial coordinate using either adiabatic potential 
energy surface (PES) obtained with Strutinsky's macroscopic-microscopic 
method [19,25] and semi-empirical approaches [23,26] or a liquid drop PES 
[20,21,22,24,27,28]. The competition between fusion and quasi-fission caused 
by the effect of the thermal fluctuactions has only been included in recent 
models of this type, e.g., [24,25,26,27,28]. The theory in [19] is completely 
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static. It only predicts the so-called cold fusion valleys in the PES that are 
formed with target-projectile combinations of closed shell nuclei. A large 
number of these combinations were succesfuUy used at GSI (Germany) to 
synthesize SHE by cold fusion reactions. The macroscopic dynamical models 
[20,21,22,24,25,26,27,28] are classical ones based on the concept of mean tra- 
jectory which is affected by the PES, the inertia and friction tensors as well 
as by thermal fluctuations. The theory of Langevin trajectories [24] is quite 
appealing, but it neglects important quantum effects (quantum fluctuations, 
ground-state shell effects, diabatic effects) and many of its ingredients are 
obtained in a rather uncertain manner, e.g., hydrodynamical inertia and 
friction tensors as well as the PES which, in the best situations [25], is cal- 
culated with a two-center shell model (TCSM) [29] that is inappropriate for 
compact nuclear shapes in very asymmetric reactions. Of course, the exper- 
imental data are not always explained. A generalized Langevin approach 
including quantum fluctuations has been proposed recently in Ref. [30]. 
• In the second type [31,32,33,34] (dinuclear system (DNS) model), the fusion 
happens only in the mass asymmetry coordinate r] = {Ai — A2)/ {Ai + A2) 
where Ai and A2 are the mass numbers of the nuclei. The DNS nuclei remain 
at the contact configuration and exchange nucleons until either all nucleons 
have been transferred from the lighter to the heavier fragment (complete 
fusion), or the DNS decays before the CN formation (quasi- fission) . The 
model assumes a sudden (double-folding in frozen density approximation) 
PES in the radial coordinate, while the PES behaves adiabatically along the 
fusion path in the 77 coordinate. Thus the theory avoids including the com- 
plex physics for large overlapping nuclei. The model is also attractive due 
to its simplicity and because its ingredients are obtained in a clearer man- 
ner (defined in terms of the properties of the individual fragments forming 
the DNS). Although it has been successfully used to explain many exper- 
imental ER cross sections, its theoretical foundation (e.g., the peculiarity 
of the PES) is not clear enough yet. Possible reasons for this could be the 
diabatic effects [35] and large values of the microscopical mass parameter 
with respect to the neck degrees of freedom [36]. 

The motivations for this work have been (i) to reconcile these conflicting the- 
ories, and (ii) to incorporate the multi-particle quantum nature of the fusing 
system, rather than assuming a continuous macroscopic fluid. The present pa- 
per is aimed at investigating the intermediate stage of CN formation and, in 
particular, the competition between fusion and quasi-flssion. In this study we 
give a quantum-statistical formulation of the problem, which can be consid- 
ered as the theoretical foundation of the exploratory calculations performed 
in Refs. [37,38]. The study is based on the following general ideas (which are 
well established but have up to now not been used in combination in any 
of the current models of fusion): (i) Once the two nuclei are at the contact 
point, the system moves in a multidimensional space of collective coordinates, 
(ii) this motion is governed by a master equation, and (iii) the nature of 
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the single-particle (sp) motion is timc-dcpcndcnt, it is initially diabatic and 
then approaches the adiabatic limit due to residual two-body collisions. As a 
result of (iii), the system moves in a time-dependent (dynamical) collective 
PES which is initially diabatic and gradually becomes adiabatic. Moreover, 
the system does not adopt a single (well-defined) shape, but as a result of 
quantum and thermal fluctuactions a probability distribution of the nuclear 
shapes [39] can be expected on the dynamical PES. The fusion is described as 
a statistically non-equilibrated evolution process of nuclear shapes, in which 
each shape evolves towards the thermal equilibrium. In addition to the dia- 
batic effects and thermal fiuctuations already treated in Ref. [37,38], we now 
include the effect of quantum fluctuations as well as rotation of the combined 
system on the shape fluctuations. The microscopic ingredients of the theory 
are obtained with a realistic TCSM based on Woods-Saxon potentials recently 
suggested in Ref. [40] . One of the merits of this TCSM is that it permits study- 
ing very mass-asymmetric reactions. The present theory is unique and more 
realistic than the current fusion theories because it is founded on microscop- 
ical grounds and also for the flrst time incorporates a wide range of physical 
effects that crucially influence the formation of the heaviest compound nuclei. 

The diabatic sp motion arises from the coherent character of the coupling 
between the intrinsic and collective degrees of freedom whilst the dissipation is 
caused by the residual two-body collisions [41]. The adiabatic limit refers to the 
case when the occupation of the sp energy levels obeys an equilibrated Fermi- 
distribution with a flnite temperature. At the diabatic limit [40,41] (elastic 
nuclear matter), the nucleons do not always occupy the lowest free sp energy 
levels as in the adiabatic case (plastic nuclear matter), but remain in the 
diabatic states during collective motion of the system. Fig. 1 illustrates this 
situation. In other words, the probability for Landau-Zcncr promotion of the 
nucleons at the avoided crossings of adiabatic sp levels is close to one. This 
approach is realistic in the approach phase of collisions near the Coulomb 
barrier where the total excitation energy per nucleon E* > 0.03 MeV [42,35]. 
This has been supported by time-dependent Hartree-Fock calculations [43]. 
During the transition from the diabatic to the adiabatic limit, the nuclear 
matter is elastoplastic like glycerine. Possible indications of the diabatic sp 
motion can be extracted from the pre-compound emission of fast particles 
[44,45,46] and gamma rays [47] as well as from the intrinsic excitation energy 
per nucleon of two quasi-flssion fragments [48] which is expected to be different 
for each fragment, i.e., the fragments would have different temperatures. 

Experimental evidence to support the point of view of a quantum-statistical 
description of the competition between fusion and quasi-flssion can be found, 
e.g., in Rcfs. [9,49]. Experiments presented in Ref. [9] show that (i) for a 
given target-projectile combination in the entrance channel many complex 
exit channels, involving large mass and charge transfers, develop after the 
passage (inside) of the Coulomb barrier, and (ii) those complex exit channels 
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adiabatic levels 



collective coordinate q 

Fig. 1. Diabatic sp motion at an avoided crossing of two molecular adiabatic sp 
levels. 

reveal that the fragments have different temperatures, i.e., the system is not 
in thermal equilibrium. The experimental study in Ref. [49] clearly indicates 
that the formation of the CN occurs in a thermally non-equilibrated process 
which depends on the entrance channel. 

In section 2 we will present the theoretical formalism in the following order: (i) 
the master equation for the evolution of the nuclear shapes is re-derived with 
both the coarse-graining procedure used by Norenberg in Ref. [50] and meth- 
ods of the non-equilibrium quantum statistical mechanics [51,52,53], (ii) the 
transition probability between the nuclear shapes is discussed, (iii) a model 
of treating the quantum and thermal effects on the shape fluctuations is pro- 
posed, and (iv) a realistic model of calculating the probability distribution of 
the nuclear shapes along with the probability for CN formation and quasi- 
fission is suggested. The calculation of other possible observables like the mass 
distribution of the quasi-fission fragments as well as the average total and 
intrinsic excitation energy of the CN is finally discussed. Here observables 
mean theoretical quantities that may be verified by direct experiments. Of 
course, the experimental data may contain contributions of many effects that 
our theory does not include yet. For instance, the experimental mass yield in 
fusion-fission reactions includes the contribution of fission events, in addition 
to what we define here as quasi-fission. We would like to emphasize that our 
theoretical mass distribution does not include the fission component (decay of 
the CN into two fragments), but is limited to all the binary fragmentations 
which occur after capture and before the CN formation. This is exactly what 
we call quasi-fission. The separation between the stage of formation of the 
CN and its fission is justified if the fission time scale is much larger than the 
formation time of the CN, which is supported by experiments in Ref. [54]. 
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In section 3, numerical results for heavy-ion reactions leading to the No 
compound nucleus are shown and discussed. Conclusions are drawn in section 
4. 



2 Theoretical formalism 

2.1 Conceptual details 

The present theory aims at describing the competition between fusion and 
quasi-fission after capture of the interacting nuclei, translated to a problem of 
diffusion of nuclear shapes through the multi-dimensional dynamical collective 
potential energy landscape. To my knowledge, there are no experimental ob- 
servations so far, which rule out this scenario. The theory is essentially based 
on the single-particle model. This may explain and predict the dependence 
of fusion, after contact of the nuclei, on the shell structure of the interacting 
nuclei. Owing to the statistical nature of this approach, there is no equation 
of motion for individual collective coordinates. The theory only describes the 
evolution in time of an ensemble of nuclear shapes (parametrically defined by 
a set of collective coordinates) that develop following contact of the interact- 
ing nuclei. The basic macroscopic variable will be the nuclear shape, which 
determines the two-center mean-field in which the nucleons are moving. After 
capture of the nuclei, the motion of the compact fusing system is expected to 
be slow (overdamped) , as the initial diabatic collective PES practically equal- 
izes the total incident energy of the system (small kinetic energy). The usual 
concept of collective inertia and forces do not appear in the theory explicitly, 
because the concept of classical trajectory is not applied. Effects similar to 
the collective inertial effects may be included in the local transition probabil- 
ity rate between the nuclear shapes. This transition probability rate clearly 
incorporates quantum and thermal effects which drive shape fluctuations. 

The master equation describing the evolution of the nuclear shapes is obtained 
using the time-dependent TCSM representation. The two-body residual inter- 
actions are responsible for the change of the occupation numbers of the two- 
center sp orbitals that define the set of many-body channel wave-functions 
(Slater determinants). The coupling between channels belonging to the same 
nuclear shape is assumed to be much stronger than the coupling between 
channels of different nuclear shapes. As a consequence, the Markov limit of 
the master equation is obtained. This is consistent with our expectation that 
the evolution in time of the compact nuclear shapes is slow. 

The phrase thermal equilibrium (local) will only refer to the situation when 
the occupation numbers of the sp orbits at a given nuclear shape obey a Fermi 
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distribution for a finite temperature. Following capture, due to diabatic effects 
in the entrance phase of the reaction, the sp occupations can be strongly de- 
viated from the Fermi distribution (there is no thermal equilibrium in the 
system and, strictly speaking, the concept of temperature is meaningless) . Di- 
abaticity only produces coherent particle-hole (ph) excitations that contribute 
to the collective PES. The two-body residual interactions gradually destroy 
these coherent ph excitations, i.e., the system heats up and the sp occupation 
numbers evolve in time towards a Fermi distribution for a finite temperature. 
The evolution in time of the nuclear shapes (a non-equilibrated macroscopic 
process) occurs in conjunction with the thermalization of each nuclear shape. 
This scenario will be modelled solving the master equation coupled to the re- 
laxation equation for the sp occupation numbers. We will not include nucleons 
evaporation. It is important to say that we cannot distinguish between coher- 
ent and incoherent ph excitations during the relaxation of the sp occupation 
numbers. For the sake of simplicity, we will assume that the initial coherent ph 
excitations turns into incoherent ones (heat) when the thermal equilibrium is 
reached, i.e., collective modes like surface vibrations of a given nuclear shape 
are completely damped. 

2.2 General aspects 

The total many-body wave-function |^(t)) of the fusing system is expanded 
in terms of a set of many-body channel functions {|$i(i))} as follows 

\^{t)) = Y.am\Ut)), (1) 

i 

which satisfies the time-dependent Schrodinger equation 

in^-^-H{t)m)). (2) 



where H{t) = T + V{t) is the total many-body Hamiltonian. The poten- 
tial V{t) contains a part which represents the two- center mean field VrcsMit) 
and the rest is related to residual nucleon-nucleon interactions K-es(^); i-e., 
V{t) — VrcsMit) + Vres{t). The evolution in time of the nuclear shapes is im- 
plicitly included by means of the time-dependent potential V{t). The channel 
functions {|$j(t))} are dynamically defined by the orthonormalization (e.g., 
with the Gram-Schmidt method) of the set of Slater determinants {|$j(t))} 
(built from the occupied two-center sp states, not necessarily the lowest or- 
bitals) of all two-center mean field configurations included, i.e., 

\Ht))-i:himM, (3) 



7 



where the unitary transformation matrix bij{t) results from the orthonormal- 
ization procedure. The coefficients aj(t) arc the occupation amphtudes of the 
individual channels which are assumed to be orthonormal 

{^j{t)\Mt))^S,j, (4) 



satisfying the relation 



Inserting (1) into (2) and projecting on the following coupled equations 

for the amplitudes are obtained 



dt .,. 



Y.a,{t)V^,{t), (6) 



where the non-diagonal matrix elements Vji{t) = h ^($-,(t)|iJ(t)|$j(t)) yield 
transitions between the many-body channels. We now define the density matrix 

Pji{t)^aj{t)a*{t), (7) 



where the diagonal elements are the time-dependent occupation probabilities 
of the individual channels. With (6) and (7), the Liouville equations for the 
evolution of Pji{t) read as 

^^PM = J2L,,,,{t)Mt), (8) 

l,k 



where Lji^ik{t) is the Liouville operator defined as 

Lji,ik{t) = Vjiit)Sik - SjiVkiit). (9) 

The operator Lji^ik{t) is a tetradic and behaves like a square matrix in the 
superspacc 7i07i*, whereas pik = {lk\p(t)) behaves like a vector in that 
superspace. 7i denotes the Hilbert space of the channel states. For details 
about the algebra of tetradics we refer to Ref. [51]. The total Hilbert space 
H is now divided into subspaces (i.e., H — ^Hn^p — 1,2,. . .N). Each 
subspace defines a macrostate with dimension rf^ which is the number 
of channels {micro states) contained in 7Y^. Each macrostate p, is considered 
to be uniquely related to a fixed geometry of the two-center mean field (i.e., 
a fixed sp spectrum of the TCSM) that is described by a set of collective 
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(shape) coordinates denoted by q (e.g., the relative distance between the two 
mean field potentials R, the so-called fragmentation coordinate 1], etc). The 
difi:erent many-body channels belonging to a given /j, are assumed to involve 
only the Slater determinants defined by those fixed sp energy levels which are 
occupied. It is important to keep in mind that each macrostate /i is assumed 
to be uniquely associated with a nuclear shape q. The occupation probability 
Pixit) of the macrostate is defined as 

P,{t) = E P^mit). (10) 
meH^ 



These macroscopic probabilities Pfj.(t) are considered to be the relevant quanti- 
ties for describing the evolution of the combined system. We shall now obtain 
an equation of motion {master equation) for these probabilities. 



2.3 Master equation 



To obtain the master equation for P^{t) we should first introduce the so-called 
coarse-graining operators [50] defined as 

~ '^/I^ E \mm){m'm'\, (11) 

m.m'eHu 



which act on the vectors of the superspace TC^ 7i^* such as the density 
matrix (7). These operators are projectors and fulfil the usual relations = 
Cn and C^Cj, = C^5^^^. The operators act on the vector \p{t)) associated 
with the density matrix taking into account the definition (10) as follows 

CMi))-d-' E \mm){m'm'\p{t))^d-'P,{t)Y.\mm). (12) 



Using the coarse-graining operators (11), the projection operators C = Y,)j.C^ 
and Q — 1 — C can be defined. With the Nakajima-Zwanzig projection tech- 
nique described in Ref. [52] and acting separately with the projectors C and 
Q on the Liouville equations (8), the following generalized master equation 
for the macroscopic probabilities Puit) is obtained (see Appendix A) 

dP it) ''f" 

^^ = E J drK,^{t,T)[d,P^{t-r)-d^P,{t-T)]+G,{t,to), (13) 

MT^*^ 



where K^^it^r) is the so-called memory kernel and G^{t,to) includes the ef- 
fects on Pu{t) due to initial correlations between the channels, i.e., non-diagonal 
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matrix elements Q\p{t(})) of the initial density matrix. The quantities Ki,^{t, r) 
(symmetric regarding the indexes v and /x) and G^it^to) read as 



K,^{t,T) = -d-' Y E {nn\L{t)eM-^Q j dt'L{t')] 

neHv meH^ 

■QL{t-T)\mm), (14) 

t 

G,{t, to) = E (^^l^(^) expHg J dt'L{t')] ■ Q\p{to)). (15) 

neHu to 



In the following, we will consider that the initial density matrix is diagonal, 

i.e. Q\p{to)) — 0, and thus the inhomogencity Gi,{t,to) in (13) is removed. It 
is certainly guaranteed if the initial time to is fixed before the decoherence in 
the process appears as will be discussed below. Hence, the generalized master 
equation (13) becomes a closed equation in the macroscopic probabilities Puit). 

Markov limit. Assuming that the coupling matrix elements Vji(t) between the 
channels are (i) essentially constant during the time interval r, and (ii) ran- 
domly Gaussian distributed [55], which is justified if many degrees of freedom 
are involved in the process as is the case in the fusion of a heavy system, the 
memory kernel (14) can be approximated as follows (see Appendix D in Ref. 
[50]) 

Tmem \t) 



where {■ ■ .)u and {■ ■ denote the mean values over the channels in the 
macrostates v and /x, respectively, and T^em(^) so-called memory time 

that determines the decay of the memory kernel. The averaging coupling and 
the memory time are defined as follows 

mUt)\%), = dZ' Y d-' E \ynm{t)\\ (17) 

neHv meHfi 

r^^iit) = v^{E((l^-(t)r). + (|V;n.(t)r).)]}-i (18) 

ieH 



The decay of the memory kernel results from the coupling of all the channels 
which destroys the coherence in the transition between the channels, i.e., the 
reversibility. The coherent character of the sp motion refers to the diabaticity of 
this motion [41]. Only during the memory time r^eml^) the transition between 
the macrostates /i and u is coherent (reversible) . It is important to emphasize 
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that the memory time is the same as the time of coherence. The coherent 
(diabatic) collective motion is expected to occur in the approach phase of 
the reaction at incident energies slightly above the Coulomb barrier [42]. The 
memory time was roughly estimated in Ref. [50] [i.e., < (8 ■ 10"^^)/A s for 
systems with mass number A, e.g., 4- 10~^^ s for A ~ 200] at energies near the 
barrier. If the macroscopic probabilities Pij,{t) are essentially constant during 
time intervals of the order of the memory times rj^^Xi'^)^ then the generalized 
master equation (13), with the memory kernel (16) and integrating over r, 
becomes the following master equation in the Markov hmit (similar to Pauli's 
master equation (see Ref. [51])) 



where Wu^it) are the so-called transition probability rate (symmetric with re- 
spect to the indexes) defined as 

^^.{t)^r^:l{t)-{{\Vnrn{t)\\),. (20) 

The master equation (19) is valid if the memory times T^em(^) represent the 
shortest time scale in the fusion process. The time scale for the local variations 
of the macroscopic probabilities P^{t) are called relaxation times r^^^i, which 
are estimated to be of the order of 10~^^ s in peripheral deep-inelastic heavy 
ion coUisions [50] at the energies considered. In addition to these time scales, 
there is a time scale associated with the recurrence of the system close to its 
initial macroscopic state that is called the Poincare recurrence time [53]. The 
Poincare recurrence time has to be the largest one for the validity of the master 
equation (19) (see Ref. [50]) and it has been estimated to be several orders 
of magnitude larger than the time scales mentioned above in deep-inelastic 
processes at the energies discussed [50] . We would hke to stress that the vahdity 
conditions of the master equation (19) in terms of the relation between the 
different time scales is less restrictive than assuming local thermal equilibrium 
in each nuclear shape. The master equation (19) seems to be justified in fusion 
of heavy nuclei near the Coulomb barrier. 

In the adiabatic TCSM the coupling between the channels can be expected to 
be weak, as in this approach the nucleons remain at the lowest free sp orbitals 
obeying an equilibrated Fermi distribution during the collective motion of the 
system. At very small temperatures, due to Pauli blocking effects, the mean 
free path of nucleons is expected to be very long. Nevertheless, the adiabatic 
two-center sp states are not used in the present work to define the many- 
body channels (3), but the time-dependent two-center sp orbits, i.e., those sp 
states that are dynamically occupied by the nucleons. Following the entrance 
phase of the collision, these orbitals are the diabatic sp states, and they turn 
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into the adiabatic ones during the fusion. My calculations with the realistic 
diabatic TCSM based on Woods-Saxon potentials [40] indicate that many 
particle-hole excitations with a considerable energy (particle states up to 10 
MeV above the Fermi surface for ^^Ca -|- ^°^Pb) can appear following capture, 
due to diabatic effects. The two-body collisions between a great number of 
energetic nucleons give rise to transitions between the many-body channels. 
The life-time of these highly excited sp states [see eq.(42)] is indeed very 
short (~ 10~^^ s) and remains so because the system heats up when the 
sp occupation numbers gradually evolve to a Fermi distribution for a finite 
temperature (1-2 MeV). These arguments supported by numerical calculations 
allow me to draw the conclusion that the coupling between the channels is 
not weak, at least within a given macrostate fi (nuclear shape), but on the 
contrary strong. As mentioned in subsection 2.1, we will assume that the 
coupling between channels belonging to different macrostates is much weaker 
than the coupling between channels within the same macrostate. As discussed 
in Ref . [50] , this strengthens the validity of the Markov approximation. At the 
moment this is the most practical approach to use in numerical calculations 
and there is no experimental evidence contradicting it. 



Saloner and Weidenmiiller quantitatively demonstrated in Ref. [56] that (i) in 
a deeply inelastic collision near the Coulomb barrier, like "^"Ar + ^^^Th at an 
incident energy of 200 MeV, the strong-coupling condition is fulfilled for the 
bulk of the interaction region, and (ii) the Markov approximation is indeed 
justified. For peripheral collisions Ayik and Norenberg discussed in general 
in Refs. [57,58] the Markov approximation in the limits of weak and strong 
coupling. In both limits the decay of the memory kernel (14) is Gaussian, 
but the memory time is different [57] . In the weak-coupling limit the memory 
time is theoretically determined by the correlation time (see eq. 4.13 in Ref. 
[57]) which may be infinite in the overdamped regime of the collective motion. 
In the strong-coupling limit the memory time (18) corresponds to the decay 
time of the many-body propagator in (14). In fusion of massive systems the 
Markov approximation (19) can be expected to be valid because the evolution 
of the nuclear shapes (where many particles participate) depends on several 
collective coordinates (deformation, mass and charge asymmetries etc) that 
may change very slowly. All slowly varying collective coordinates are included 
in the definition of the nuclear shapes that evolve obeying the master equation 
(19). To summarize, we believe that after a short period of coherent motion 
(memory time or time of coherence) during the approach phase of the reaction, 
the system gradually evolves to a local thermal equilibrium and the probability 
distribution of the nuclear shapes may be described by the master equation 
(19). Now we will discuss the transition probability rate between the nuclear 
shapes. 
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2-4 Transition probability rate 



Inserting expressions (17) and (18) into expression (20), the transition prob- 
ability rate w^i^{t) can be re- written as 



where 



X7{t)^V2^- Y: E \Vnm(t)\'- 
neHv meT-Lfi 

{EK- E iKi(^)r+4- E iK^iWri}"'- (22) 

ieH neHv meHu 



Looking at the master equation (19) we can see that net transition probabihty 
rate K^^it) (not symmetric regarding the indexes) between the macrostates 
and u can be defined as follows 

K,{t) = w.,{t)-d,^\7{t)-[^]K (23) 
K,At) = w,,{t)-d, = \^''{t)-[^]^. (24) 



We would hke to note that the transition probability rate (21) and the net tran- 
sition probabihty rates (23) and (24) adopt similar forms as those proposed 
by Morctto and Sventek in Ref. [59] based on phenomenological arguments 
related to the level density of the system. This will be demonstrated in the 
remainder of this subsection. 

In the phenomenological expressions in Ref. [59] the quantity Ao^(t) is a nuclear 
shape-dependent geometrical parameter and the dimension of the macrostates 
is replaced by their level densities q. The dimension of the macrostates {d^ 
or dfj) can be expressed in terms of their entropies following the definition of 
entropy of the statistical mechanics, i.e. 3,^ ~ \n.{di,). Thus d^, ~ exp(5',^) ~ 
Q{E*{t)), being E*(t) the total excitation energy of the macrostate u. Hence, 
the net transition probability rate, e.g. A,^^(i), can be written as 

A..(^) = Ao^'^(i)-[^|^]^ (25) 
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where the proportionahty coefficient between the dimension of the macrostates 
and the level density is assumed to be a constant independent of the macrostate 
{di, = AEg{E*{t)), where AE = h/Trei being Tj-ei an average of the relaxation 
times r^g; associated with the macroscopic probabilities Pu{t) discussed above). 
The first factor in (25), Xo^{t), essentially includes the quantum effects on the 
evolution of the macrostates {shape fluctuations) due to the coupling between 
the channels, whereas the second factor is solely related to thermal cSccts. The 
quantum effects arc actually not completely separated from the thermal ones 
in expression (25) owing to the dependence of Ao^(t) on the dimension of the 
macrostates [see expression (22)]. If the Fermi gas formula for level density is 
used, i.e., g{x) ~ exp{2^/ax) being a — A/12 MeV"^ (^4 is the total mass of 
the system), then expression (25) reads 

K,{t) = K^{t) ■ exp[^aE*it) - sJaE;{t)]. (26) 

In the limit of a thermally equilibrated collective motion, the available total 
excitation energy (shared by the intrinsic and collective degrees of freedom) 
can be expressed as El{t) — E — V^^^^^f,, being E the total energy (measured 
from the deepest minimum in the potential energy) and V^'^j^jj the so-called 
adiabatic potential energy. Performing a Taylor expansion of In f?(-E — V^'^j^j) in 
powers of Vj^j^fe and taking contributions up to the first term, the level density 
g{E*{t)) fti g{E) exp(— V^^^^j^^/T), where T is the nuclear temperature defined 
as ~ dln[g{x)]/dx\x=E [59]. Inserting this result into expression (25), the 
well-known expression of the net transition probability proposed by Moretto 
and Sventek in Ref. [59] is obtained, i.e., A^^(t) ~ ^w[(y^diab ~ ^adiab)/'^'^]- 
Most of the applications so far [26,38,59,60] have used this approach in which 
only the thermal effects on the shape fiuctuactions were included. 

It is worth mentioning that a semi-microscopical description of the transi- 
tion probabilities exits [61,62], but it is valid only if the overlap between the 
nuclei is small. However, the collective motion, which we are considering, is 
neither thermally equilibrated nor peripheral and, therefore, the net transi- 
tion probabilities may have to be calculated in a different way. In the next 
two subsections, we show a model to calculate the factors Xg'^lt) and the total 
excitation energy E*{t), respectively. 

2.5 Model for calculating the factors Ao'*(i) 

To calculate Xo^{t) in expression (22), a modelling of the non-diagonal couphng 
between the many-body channels, e.g., Vnm{t), is needed. Using the diabatic 
two-center sp basis [40] (which minimizes the radial couplings) and neglecting 
the rotational (or Coriolis) couplings (which may be justified if the system 
is heavy), this coupling matrix element is essentially given by the coupling 
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between the Slater determinants caused by the residual interaction 

Vres{t), i.e., 

^ E brnfim^m^AWresmAt)), (27) 

./ ./ 

« J 

provided the transformation matrixes b^j'{t) and b^^'{t) associated with the 
macrostates fi and u, respectively, are known (the channel functions |$,„(t)) 
and |$„(t)) as well as the set of Slater determinants {|$j'(^))} and {|$j'(t))} 
belong to the macrostates /x and respectively). 

We now assume that the two-body residual interaction only depends on the 
relative position between the micleons, i.e., K-es(^) = 'I2i>jv{ri — Vj). The 
matrix elements among the Slater determinants in expression (27) will be 
denoted by V-'j' (t) and they are calculated as follows 

i>j 

(^a,_i(ri-i)|^ft_i(ri_i))... 
(<^a,_i(rj-i)|<?/3,_i(ri-i)) • 
{<Pai{ri)<Paj{rj)\v{ri-rj)\{pf}.{ri){pp.{rj)) ■ 

(</?ai+i(ri+i)lm+i(ri+i)) ... 

(</'a,+i(rj+l)|<^^,.+i(r^-+l)) . . . {iPaA^A)\^f3A^A)) }, (28) 

where ipa and ipjs are the occupied diabatic molecular sp states with the set of 
quantum numbers a and (3, which build up the Slater determinants (t)) and 
\^j'{t)), respectively. A'^a) ^^'^ '^(/S) refer to the antisymmetrization operators 
in the macrostates u and /i, respectively, which act on the set of quantum 
numbers a and P of the sp states. These operators are defined as follows, 
e.g., .A^'q.) = I]p(— l)'^*^^^ . . ., where P denotes a permutation of the quantum 
numbers a among the nucleons (i.e., P{1) — cki, . . . P{A) — a^) and 7r(P) 
corresponds to the number of inversions (even or odd) of the permutation P, 
i.e., (-l)-(^) = ±1. 

Using in (28) the antisymmetrization operator A'^f^^ and considering the dia- 
batic feature of the sp basis, only diagonal terms appear in the overlap matrix 
elements of the sp states. Demonstrated in Fig.2 using the TCSM [40] and in 
Ref. [55] with the Nilsson-model, the average of these diagonal overlap ma- 
trix elements can be fitted with a local Gaussian with a variance a^, e.g., 
(<yi'ai(ri)|<^ai(ri)) ^ exp(-(q^ - q^f/2al), where and are the collective 
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Fig. 2. (Color online) Diagonal overlap of two-center sp states with the projection 
= 1/2 of the total angular momentum for the symmetric fragmentation of ^^^No 
around R = 12 fm. The thick solid curve shows a Gaussian fitted by the method of 
least squares. 

(shape) coordinates associated with the macrostates z/ and /i, respectively. 
Please note that for reasons of clarity we refer to a Gaussian with a vari- 
ance, but it may be a product of Gaussians with different variance where each 
of them is related to an individual collective coordinate in the set q. Thus 
expression (28) becomes 



II^M^ {Vo.Ari)'^a,{rj)\v{ri - rj)|<^a,(ri)<^a^.(rj))}. (29) 

i>j 

Assuming that f (r^ — r^) = 5(rj — r^) (zero-range interaction), then 

{y^aA^i)<^a,{rj)\v{ri - rj)\^^X^i)(p^^{rj)) = (y?«^(ri)y?«^(ri)|^^^(ri)(^^^.(ri)) 

= Qi- (30) 
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Denoting by Clf the antisymmetrization of (30) in (29) including the factor 



{A\) ^, the matrix elements (29) are finally written as 



^2 , , 



y,,.(t) ^ exp{-i^L^} . Y^cl/ , (31) 

i>j 



where a^, — a^/ \J A — 2. Thus the matrix elements (27) are expressed as follows 

Km(i) ~ ^^V{-^^^^^^} ■ Cnm. (32) 



where 



Cnm = Y.bmf{t)KAi) (^Clf)- (33) 



-■' -■' i>j 



If we now use expression (32) in expression (22), the quantities Xo^{t) can be 
written as follows 

X7{t) « Ko'it) ■ exp{-^^^LS^l)!}, (34) 

(Tf, 



where 



neHv meH^ 



leH neHv meHu 



(35) 



The calculation of the dimensions dy and in terms of the level density was 
discussed in the previous subsection. 



2.6 Model for calculating the total excitation energy E*{t) 



The available total excitation energy E*{t) > in the rotating body- fixed ref- 
erence frame, which can be distributed between the intrinsic and the remaining 
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collective degrees of freedom, is calculated as follows 

Etit) = Maximum {e^i), E,,m. - ^^'^^/q^ - Vdyn{t)}, (36) 

where e*(i) is the intrinsic excitation energy, E(.,m. is the total incident energy 
in the center of mass reference frame, J is the total angular momentum of 
the system with a local moment of inertia Q^, and VJy^{t) is the dynamical 
collective PES defined as 

V^ynit) = Kd^a6«(i)) + ^V^iM, (37) 

where K'^jab is the adiabatic PES which is the sum of the liquid drop energy 
and the microscopic shell and pairing corrections obtained with Strutinsky's 
method. The so-called diabatic contribution AVJ^^^(t) is initially maximal, but 
gradually decreases when the sp occupation numbers approach the equili- 
brated Fermi distribution. The dynamical PES (37) describes a continuous 
transition from the initial diabatic potential to the asymptotic adiabatic one. 
We would hke to stress that the effect of the rotation on the evolution of the 
combined system will only be included through expression (36). This expres- 
sion clearly shows that in a classically inaccessible macrostate u [i.e., the sec- 
ond part of (36) would be negative] the intrinsic excitation energy determines 
E*{t). In the remaining part of this subsection, we discuss the calculation of 
the diabatic contribution Ayj<g^(i) and the intrinsic excitation energy e*(i). 

Diabatic contribution. According to the dissipative diabatic dynamics (DDD) 
approach (see [47] and references therein), which is consistent with the fea- 
tures of the collective motion we aim to describe, the dissipation arises from 
the decay of the diabatic contribution AV^-^^f^{t) due to the residual two-body 
interactions (neglecting the effect of the rotational couplings). The diabatic 
contribution AV'J^^b(t) results from the incident collective kinetic energy pri- 
marily stored in the system as a reversible potential energy (elastic response) 
during the initial coherent phase of the sp motion. Following this coherent 
stage, this amount of energy gradually turns into heat or intrinsic excitation 
energy, i.e., a thermalization process occurs. Using the adiabatic and diabatic 
sp energy levels of the TCSM [40,63], AV'J^„b(t) can be calculated as follows 
[63] 

^V^ia^it) = E • K^'^'it) - E • nf'^\e^, n (38) 

" a 

where a and d denote the quantum numbers of the diabatic and adiabatic 
levels, respectively. The diabatic contribution of neutrons and protons is sepa- 
rately calculated. The adiabatic quantum numbers are contained in the set of 



18 



the diabatic ones [64]. The diabatic sp occupations are indicated by 
while the adiabatic occupations by n|^f'^"*^(ei?, T). The adiabatic occupations 
obey a Fermi distribution for a chemical potential ep{t) and a temperature 
parameter T{t) that only measures the local intrinsic excitation energy e*(t). 
Please note that the time dependence of the chemical potential and the tem- 
perature parameter is not included in the expressions of this subsection for 
reasons of clarity. For a given temperature parameter the chemical potential is 
obtained by the conservation of the particle number. The calculation of T{t) 
will be explained below. Since the diabatic levels solely differ from the adi- 
abatic levels at some of the avoided crossings of the adiabatic states [40,63], 
an approximation to the time-dependent diabatic contribution (38) can be 
written as [35] 

AV^L^l^) ^ E ■ K""'\t) - nT'-^'icF, T)]. (39) 



The time-dependent diabatic occupation numbers n[^''^*"^(t) (actual occupa- 
tions within this model) are described with the following relaxation equation 

[42] 

= -^."'(^) • K'^"'(^) - <'"'^"'(e^, T)l (40) 



where is an average relaxation time (in order to conserve the number of 
particles). The initial diabatic occupations nJ^'*"''(to) are the same for all 
the macrostates (shapes) which arise from some equilibrium distribution 
n^°'"*"*(To). Here vq denotes a shape of the entrance system before the dia- 
baticity is manifested, e.g., the nuclei are well-separated at their ground-states. 
The parameter Tq will only be used to simulate a diffused distribution of sp oc- 
cupations around the Fermi surface of the individual nuclei, i.e., the step func- 
tion is somewhat smeared out. This can happen although the nuclei are cold 
due to the pairing correlations. Since initially in the entrance channel we deal 
with two separated nuclei. To is in general different for each nucleus. There- 
fore, To will be denoted as a two-component parameter, i.e.. To = (Toi,To2). 
Realistic values of To can be obtained solving the gap equation [65] for each 
isolated nucleus. The relaxation time Ty is defined as 

1 _ Ea<^°''(t)- r-(e^,T) 



where the widths FJ^ of the sp levels due to residual interactions are obtained 
with the parametrization given in Ref. [66] 

K = To ^(6^''^^"" - epf + (7rT)V(l + [{e''/"'' - e^f + {iiTfyc'). (42) 
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where O.OSMeV"^ ^ Tq ^ ^ O.OGlMeV-^ (depending on the nuclear density) 
and 15MeV^ c ^ 30MeV. 

Using expression (39) and equation (40), the diabatic contribution A\/j<^^(t) 
satisfies the following differential equation 

^^^^ + :^^^VL,{t) = 0, (43) 



whose solution reads as 

t 

^VLbit) = ^VLM ■ exp[- / dt'T-\t')]. (44) 

to 



Intrinsic excitation. The temperature parameter T{t) is related to the local 
intrinsic excitation e*(t) through the Fermi gas expression e*{t) — a ■ T{tY, 
being a — A/V2 MeV~^ (^4 is the total mass of the system). For the calculation 
of the local intrinsic excitation energy, we follow the dissipation mechanism of 
the DDD approach and make the following ansatz 



if Ec.m. — ^ 29"*"^^ ~ ^dyni^) — 0) othcrwisc the intrinsic excitation e* (t) remains 
constant. The system only heats up if it is in a classical region. The solution 
of (45) taking into account expression (44) is as follows 



t tl 
el{t) = e:(to) + ^V^iaM ■ I dtir-\ti) ■ exp[- J dt2T-\t2)]. 

to to 



(46) 



The solution (46) is valid following the coherent period of the reaction. The 
initial time to will therefore be fixed just after the coherence period (memory 
time Tmem), whcu the master equation (19) also applies. The initial diabatic 
contribution AVj^g^(to), calculated with expression (39), arises from the pre- 
vious period of coherence (diabaticity) of the sp motion which started with 
the well-separated nuclei having some equilibrated Fermi distribution for the 
sp occupations. Since during the coherence period the system does not warm 
up, the initial intrinsic excitation energy of the macrostates e*(to) is equal to 
^loito — Tmem) ■ We wiU assumc that the well-separated colliding nuclei are cold 
[C^ito - Tmem) ^ MeY]. 
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Since at the initial time to some nuclear shapes (members of the ensemble) 
may be in a classically forbidden region, they start heating up later on, at the 

time to > ^Oj when the condition i?c.m. — ^ 20'^^'' ^ ^dyni^o) = is fulfilled due 
to the relaxation of the initial diabatic PES. The initial diabatic contribution 
Ayj<„j,(to), used in expression (46), is then calculated with expression (44). 

2.7 Model for calculating the probability distribution of the nuclear shapes 

The master equation (19), the net transition probability rate (26) along with 
the Ao^(t) factors (34) and the total excitation energies (36) build up a re- 
alistic model to study the probability distribution of the nuclear shapes dur- 
ing the fusion process, i.e., the calculation of the macroscopic probabilities 
Pi,{t). First we will describe the calculation of the probability for CN for- 
mation PcN and the quasi-fission probability Pqf, which result from the 
time-dependent probability distribution of the nuclear shapes Pvif). After- 
wards we will discuss the calculation of other possible observables like the 
mass distribution of the quasi-fission fragments as well as the average total 
and intrinsic excitation energy of the CN. It is important to keep in mind that 
all these quantities are referred to the rotating body-fixed reference frame, i.e., 
the fusing system has a constant total angular momentum J. 

2.8 The probability for CN formation and quasi-fission, the mass distribution 
of the quasi-fission fragments as well as the average total and intrinsic 
excitation energy of the CN 

The space of the macrostates or nuclear shapes is divided into three re- 
gions: (i) region of compact shapes around the spherical shape of the CN 
{fusion region), (ii) region of separated fragments beyond the Coulomb bar- 
rier [quasi-fission region), and (iii) region of intermediate shapes which could 
lead to fusion or quasi-fission {competition region). For each fragmentation 
T] — {Ai — A2)/{Ai -\- A2), being Ai and A2 the mass of the nuclei, the fusion 
radius Rp is defined as the distance from which the shell structure of the spher- 
ical CN manifests. Studying the TCSM levels diagram [40] we found that at 
the radius Rp the lowest sp energy level (with Vt = 1/2 being the projection 
of the sp total angular momentum on the internuclear axis) reaches a minimal 
value, e.g., see Fig. 3 where the adiabatic neutron (bottom) and proton (top) 
levels correlation diagrams for the reaction ^^Al -|- ^^^Ac — > ^^^No (77 = 0.75) 
are shown. In this example, Rp is about 6 fm. Fig. 4 shows the fusion radius 
Rp (black solid curve), the contact radius Rc (blue dotted curve) and the 
radius of the adiabatic Coulomb barrier Rb for a central collision (red solid 
curve) as a function of the entrance channel 77 leading to ^^^No. The radius 
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Fig. 3. (Color online) The lowest adiabatic neutron (bottom) and proton (top) 
molecular sp levels as a function of the radius R between the nuclei for the reac- 



tion ^^Al + Ac 



256 



No. Different curves correspond to different values of the 



projection 0, of the sp total angular momentum on the internuclear axis. See sect. 
3 for further details. 

Rb may change as a function of the total angular momentum J. In sect. 3 
further details are given. 

Initially it is assumed that the system is at the contact configuration (in the 
competition region inside the Coulomb barrier) with probability 1. During 
the evolution of the combined system, the probability distribution of the nu- 
clear shapes spreads over all the regions and, therefore, the total occupation 
probability for each region can be calculated as the sum of the macroscopic 
probabilities Pu{t) belonging to each region. The time- dependent probability 
for CN formation and quasi-fission are defined as the sum of the macroscopic 
probabilities Pu(t) in the fusion and quasi-fission regions, respectively. These 
read as 



ueFus. 

PQF{t) = E ^^W- 
ueQf. 



(47) 
(48) 



The time scale r^/ for fusion and quasi-fission is obtained (with a given accu- 
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Fig. 4. (Color online) Radius between the nuclei for an entrance channel t] leading 
to ^^^No, which defines the fusion {Rp, black solid curve) and the quasi- fission {Rb, 
red solid curve) regions. The dotted curve in the competition region denotes the 
contact radius Rc- See text and sect. 3 for further details. 

racy) from the condition 



PcNiTgf) + PQF{Tgf) ^ 1, 



(49) 



i.e., finally the probability for nuclear shapes in the competition region van- 
ishes. PcNijqf) and PQpijqf) are the fusion and quasi-fission probabilities, 
respectively, observed in the reaction. The observed mass distribution of the 
quasi-fission fragments Pqpijq^Tqf) can be calculated using Pyijqf) of the nu- 
clear shapes belonging to the quasi-fission region as follows 

PQF{ri,Tqf)= Y: P.iv,rqf), (50) 

ueQf. 



Please note that u in expression (50) refers to collective (shape) coordinates 
other than the fixed mass asymmetry t], e.g., the radius R between the frag- 
ments. Finally, the average total excitation energy of the CN, i.e. EQj^{Tqf), 
can be defined as 

7-T* \ YlutFus. Pvi^qf) ■ Ejj{Tqf) 

^CNK^qf) = p / N • (51) 
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In a similar way, the average intrinsic excitation energy of the CN, i.e. ^cN{'^qf)i 
can also be calculated with an expression like (51) just replacing E*{Tqf) by 



3 Numerical results 

3. 1 General features and methods 

The basic ingredients of the model are the sp spectra for the different nuclear 
shapes and the initial occupation number of the sp levels. The thermal effect 
on the sp energies can be neglected at low temperatures (T < 5 — 6 MeV) 
[67], therefore the sp levels are those of a cold system. The sp spectra for 
the different nuclear shapes are calculated solving the two-center problem [40] 
for neutrons and protons with realistic Woods-Saxon (WS) potentials, while 
the initial sp occupations are obtained with the diabatic sp motion in the 
entrance phase of the collision. The WS potential of the separated spherical 
nuclei and the spherical CN are obtained with the global parametrisation given 
by Soloviev in Ref. [68]. In a realistic calculation, the potential parameters 
should be selected in such a way that the experimental sp energies around 
the Fermi level are reproduced. The two-center potential of the fusing system 
is calculated applying the condition of volume conservation [40]. In order to 
reduce the computation time in calculating the sp spectra, a small number 
of harmonic oscillator basis states [Imax — 7 (the number of partial waves in 
which the potential acts) and n,nax = 1 (the number of separable terms in each 
partial wave)] will be included. As collective coordinates we use the radius R 
between the two WS potentials (describing the radial motion of the nuclei) 
and the asymptotic mass asymmetry coordinate rj (describing the different 
fragmentations). We will assume mass/charge equilibrium in the system, i.e., 
the charge asymmetry coordinate r]z — {Zi — Z-2)I{Z\ -|- Z-i) is equal to 77. 

Fig. 5 shows nuclear shapes of the system ^^^No as a function of R and r]. 
The shapes with negative r\ values are the same, but the position of the nuclei 
is reversed. We will use a positive value of t] for the entrance channel. From 
i? = up to the contact radius Rc = 1.2 ■ (A]^'''^ + Ag'''^) fm of the two spherical 
fragments = ^^(1 + rj) and A2 = ^(1 — rj), the forms are assumed to be 
the same as the equipotential shapes defined by the neutron Fermi surface of 
the fused system. For the sake of simplicity, we have assumed that the nuclear 
shapes for neutrons coincide with those for protons. The volume of all these 
shapes is equal to the volume of the spherical ^^^No compound nucleus. It is 
important to notice that our family of nuclear shapes in Fig. 5 does not include 
shapes with an elongated neck near the contact radius, which are commonly 
used for describing the fission process. This may not significantly affect our 
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Fig. 5. Nuclear shapes of the system ^^^No as a function of R and r]. See text for 
further details. 

results if the fission valley is mainly reached through compact fused shapes. 
In a recent work by Ichikawa et al. in Ref. [69] it was shown for cold fusion 
reactions that a ridge of several MeV, around the contact radius, separates 
the fusionlike valley from the fission valley. The height and persistence of that 
ridge for compact shapes may influence the CN formation. Of course, this 
issue needs further investigation. As mentioned in the introduction, we do not 
include fusion-fission events in the present calculations. Moreover, we will not 
include quasi-fission events that may happen through the fission valley. 

Using the axial symmetric shapes like those in Fig. 5, the nuclear part of the 
liquid drop energy (LDE) in the dynamical collective PES (37) is calculated 
with the Yukawa-plus-exponcntial method [70]. For the separated nuclei the 
LDE is equal to the Krappc-Nix-Sierk (KNS) potential [70]. The KNS poten- 
tial very well agrees with the recent semi-microscopical potential calculated by 
Denisov and Norenberg in Ref. [71] using the formalism of the energy-density 
functional along with an extended Thomas- Fermi approximation. The shell 
corrections are obtained with the method recently proposed by the author 
in Ref. [72], while the pairing corrections to the PES will be neglected. The 
dependence of the shell corrections on the intrinsic excitation energy e*(t) is 
included as usual, i.e., by using an exponential damping factor exp(— 76*) be- 
ing 7 = (1 + 1.3A-^/^)/{5A8A^/^) MeV-^ [73]. The moment of inertia of the 
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compact nuclear shapes for the rotation around an axis that is perpendicular 
to the axial symmetric axis is calculated using the rigid-rotor approximation. 
For the separated fragments, the moment of inertia is that corresponding to 
the orbital motion, i.e., ixR^ being ii the reduced mass for a fragmentation 77. 
We will study the dependence of the observables on the angular momentum 
J of the combined system. 

Since the diabatic sp excitations occur around the Fermi surface, the values 
Fq ^ = 0.061 MeV~^ for half saturation density and c = 20 MeV can be used 
in expression (42) as standard values for calculating the width of the sp lev- 
els along with the local relaxation times (41). It can be anticipated that the 
calculation very weakly depends on the parameter c, while the dependence on 
Fq ^ will be studied below. The dependence of the calculation on small values 
of To simulating the pairing correlations in the entrance system will not be in- 
vestigated in the present work, therefore a step function [Tq = (0.0, 0.0) MeV] 
will be adopted for the initial sp occupations around the Fermi surface of the 
colliding individual nuclei. Besides these parameters involved in calculating 
the total excitation energies E*{t), the first term of the XQ^{t) factors (34), i.e. 
Ko'^(t), will be regarded as a constant parameter hq that is independent of the 
macrostates. This approach means that only the gross feature of the quantum 
effects on the shape fluctuations will be included. Experiments related to the 
quasi-fission of heavy systems with Z = 108, which were analyzed [60] with 
a diffusion model based on both the master equation (19) and net transition 
probabilities like those by Moretto and Sventek discussed in subsect. 2.3, in- 
dicate that a realistic value of Kq is of the order of 10^^ s~^ hke the magnitude 
of the factor v^^"^ in (35). We will study the dependence of the calculation 
on kq values around this standard value. 

The local variance of the Gaussians in (34), i.e. and a^, can be obtained 
by means of the overlap of the TCSM sp wave-functions as shown in Fig. 2. 
However, in order to simplify the calculation, these variances will be considered 
as constants (5/? and cr^J that only depend on the region (fusion, competition 
and quasi-fission) of the nuclear shapes. Moreover, we will assume that is 
the same in all regions, whereas adopts the same value for overlapping nuclei 
(fusion and competition regions) and it is essentially zero in the quasi-fission 
region. For well-separated nuclei, as the shapes belonging to the quasi-fission 
region, the motion in the mass asymmetry coordinate 77 (nucleon transfer) is 
expected to be negligible. Unless other values of and (Jr arc pointed out, 
we will use as standard values in the calculations = (0.2, 0.2, 10^^) (each 
component refers to the fusion, the competition and the quasi-fission region, 
respectively) and aji — 0.5 fm. 

The set of equations (19), (26), (34) and (36) are solved by successive iterations 
using a small time step At = 10^^^ s. The master equation (19) is solved on a 
2D-mesh (Fig. 6) defined by the coordinates R and 77 (i? < 16 fm and \ri\ < 1 
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Fig. 6. Schematic illustration of the 2D-mcsh used in the calculation. A node of 
the mesh is related to a nuclear shape. The thick solid curves separate the fusion, 
competition and quasi-fission regions. See text for further details 



with steps Ai? and A77 defined below) . The step values Ai? and A?7 should be 
smaller than the Gaussian variances and a^, respectively, so that the mesh 
points can be interconnected. As standard values in the calculations below 
the values Ai? = 0.27 fm and Ar] = 0.1 will be used. At = ±1 (spherical 
CN) only the mesh point with i? = is included because only this point can 
be physically realized. For the sake of simplicity, as initial condition we will 
assume that Pvito) — S^uo denoting z/q the contact configuration of the entrance 
system, i.e., the system is well located at the contact configuration. Hence, we 
will neglect shape fluctuations that might happen during the capture stage. 
This effect may be included using a Gaussian distribution function for Pvito) 
instead of a delta function. The boundary conditions are defined in such a 
way that the initial Pi/ (to) irreversibly flows to the fusion and the quasi- flssion 
regions. All mesh points (shapes) belonging to the fusion and the competition 
region, respectively, are interconnected. In the quasi-flssion region, where the 
nuclei are well-separated, only the decay of the system along R occurs. The 
accuracy for expression (49) will be about 99% . 
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Fig. 7. (Color online) Entrance diabatic PES as a function of R and rj for different 
entrance channel mass asymmetries % (arrow) leading to the system ^^^No. See 
text for further details. 

3.2 Performance of the model 



3.2.1 Entrance diabatic and asymptotic adiabatic PES 

Fig. 7 shows the entrance diabatic PES as a function of R and rj for different 
mass asymmetries r)o of the entrance channel leading to the ^^^No compound 
nucleus. The energy is normalized with respect to the LDE of the spherical 
^^^No. Fig. 8 (top) shows the adiabatic PES of the cold system (the ground- 
state shell corrections are included), while at the bottom in Fig. 8 one can 
observe the LDE only, where the quantum shell effects are removed. The warm 
fusing system may reach an asymptotic adiabatic PES between these two PES 
after a gradual relaxation of the entrance diabatic PES of Fig. 7 during the 
reaction. Fig. 9 shows the diabatic capture barrier (black solid curve) for 
the entrance channels shown in Fig. 7 along with their reference adiabatic 
barriers including shell corrections (red solid curve) as well as the LDE only 
(blue dashed curve). The arrow indicates the contact radius Rc- In Fig. 10 
we show the so-called driving potential (defined as the potential energy of the 
dinuclear system as a function of r] at the contact radius Rc) arising from the 
entrance diabatic PES of Fig. 7 (curves other than the two lowest ones which 
are also presented in the small figure inserted, more details appear in Fig. 10) 
and from the asymptotic adiabatic PES of Fig. 8 [red solid curve (with shell 
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Fig. 8. (Color online) Adiabatic PES as a function of R and r/ for the cold system 
256]\Jq including ground-state shell corrections (top) and the macroscopic LDE only 
(bottom). The warm fusing system may reach an asymptotic adiabatic PES between 
these two during the reaction. See text for further details. 



corrections) and black solid curve (LDE only)]. The difference between the 
entrance diabatic driving potential and the asymptotic adiabatic one is due 
to the diabatic contribution ^Vdiab (37). From studying Figs. 7-10, we draw 
the following conclusions: 

• The entrance diabatic PES show a repulsive core at small radii R as well 
as at large mass asymmetries rj, which decreases with increasing entrance 
channel asymmetry 779. 

• The diabatic effects generally cause a capture valley in the PES, which is 
more pronounced with increasing asymmetry r^o of the reaction. 

• The entrance diabatic driving potential largely differs from the asymptotic 
adiabatic one, so the entrance diabatic PES does not completely justify the 
PES assumed in the DNS-model [31,32,33,34] of fusion. 

• The entrance diabatic PES show structures caused by the diabatic sp mo- 
tion. It is not related to the well-known (static) ground-state shell cor- 
rections to the macroscopic LDE, but to the sp motion through the shell 
structure of the different nuclear shapes. 

• First shown by Berdichesky et al. in Ref. [74], the diabatic effects raise the 
capture barrier (diabatic shift) and slightly reduce the barrier radius, which 
depends on the entrance channel. This diabatic shift seems to decrease with 



29 




10 11 12 13 14 15 16 10 11 12 13 14 15 16 



R (fm) R (fm) 

Fig. 9. (Color online) Entrance diabatic capture barrier (black solid curve) for the 
reactions shown in Fig. 7 along with their reference adiabatic barriers (red solid 
curve) and the LDE only (blue dashed curve). The arrow indicates the contact 
radius. See text for further details. 

increasing entrance channel asymmetry r]o. In some reactions a diabatic 
screening of the adiabatic potential pocket can occur. Indeed the present 
results confirm that the diabatic effects in the capture stage of the reaction 
should play an important role in the onset of fusion hindrance for heavy 
systems. This will be thoroughly studied in a separate paper. 



3.2.2 Probability distribution of the nuclear shapes and their intrinsic exci- 
tation energy 

Figs. 11-12 show the probability distribution of the nuclear shapes of the 
system ^^^No as a function of time for an entrance channel mass asymmetry 
rjo = 0.0 and 0.625, respectively. For this calculation standard values of the 
parameters have been used. The total incident energy and the total angular 
momentum are -E'c.m. = 30 MeV and J = h, respectively. It can be observed 
that initially during a certain period of time (few units of 10~^^ s) the maximal 
probability remains around the contact configuration, spreading slightly in the 
direction of symmetric fragmentations and later on, it becomes well spread out 
over compact shapes in the fusion region. The initial diabatic PES practically 
relaxes into the adiabatic PES in a time about 5*10~^^ s, and before this period 
of time, mainly dinuclear configurations near their contact point develop and 
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Fig. 10. (Color online) Driving potential resulted (i) from the entrance diabatic PES 
of Fig. 7 (curves other than the two lowest ones which are also shown in the small 
picture inserted), (ii) from the adiabatic PES at the top of Fig. 8 (red solid curve), 
and (iii) from the liquid drop PES at the bottom of Fig. 8 (black solid curve). See 
text for further details. 

decay. Afterwards, when the repulsive core of the PES at small R no longer 
exists, fused compact configurations are more probable. The shell corrections 
to the adiabatic PES (discussed below) remain and play an important role in 
the final evolution of the nuclear shapes after the dissolution of the diabatic 
contribution to the PES. This can be seen in the structures observed in the 
fusion region in Figs. 11-12 (see also Fig 8). 

Figs. 13-14 show the evolution of the intrinsic excitation energy of the nu- 
clear shapes in Figs. 11-12. Very black regions are classically forbidden. Nu- 
clear shapes there start heating up when these regions become classically al- 
lowed. At the quasi-fission time Tgf = 2 * 10~^^ s the compact fused con- 
figurations are notably warmer than the quasi-fission configurations. There- 
fore, the remaining ground-state shell effects on the PES are larger for the 
quasi-fission configurations than for the compact fused configurations. Frag- 
mentations around the entrance channel mass asymmetry rjo remain quite 
cold (quasi-elastic events), while the intrinsic excitation energy of the frag- 
mentations far away from rjo increases. However, some maxima can clearly be 
observed such as those around r] ~ ±0.4 [€cNi'''qf) ~ 25 MeV] for r]o = 0.0 
(Fig. 13) and rj ~ 0.0 [^cNi^qf) ~ 30 MeV] for rjo = 0.625. These structures as 
well as those observed in the fusion region are due to the diabatic effects and 
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Fig. 11. (Color online) Probability distribution of the nuclear shapes as a function 
of time for an entrance channel mass asymmetry r]Q = 0.0 leading to ^^^No. See text 
for further details. 

the shell corrections in the dynamical PES (see Figs. 8 and 10). Although the 
fusion-fission events are not included in the present calculations, it is shown 
that the quasi-fission fragments in the symmetric mass region for the entrance 
channel rjo = 0.625 (very close to the system ^'^Ca + ^°^Pb) are notably much 
warmer (about a factor two) than those quasi-fission products with t] ~ ±0.4 
(around the fragmentation ''°Ni + ^^^W) that correspond to a local maximum 
in the quasi-fission mass distribution discussed below. This is consistent with 
the experimental data of neutron and gamma multiplicity of the quasi-fission 
and fission products for the reaction ^^Ca + ^^'^Pb discussed by Itkis et al. in 
Ref. [15]. 

In the fusion region, the intrinsic excitation energy of the nuclear shapes clearly 
correlates with the distribution of probability shown in Figs. 11-12. For com- 
pact fused shapes, the intrinsic excitation energy coincides with the total ex- 
citation energy. Just after fusion, the average intrinsic excitation energy of 
the CN e*cj^{Tqf) is about 36.72 MeV for r]o = 0.0 in Fig. 13 and about 36.76 
MeV for rjQ = 0.625 in Fig. 14. Please note that the total incident energy -Ec.m. 
like the PES is normalized with respect to the LDE of the spherical ^^^No, so 
^CNi^qf) > Ec.m. due to negative values of the microscopic shell corrections. 
The permanence of the shell corrections in the PES is crucial for the survival 
of the heaviest CN against fission. The formation of ER will be analysed in a 
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Fig. 12. (Color online) The same as Fig. 11, but for r]Q = 0.625. See text for further 
details. 

future work. 

3.2.3 Dependence of the observables on the parameters 

We will now concentrate on (i) the convergence of the observables with respect 
to the density of the mesh (i.e., AR and Arj) for fixed values of ctr and a^, 
(ii) the dependence of the calculation on the model critical parameters (Fg ^ 
and kq), and (iii) the physical behaviour of the observables regarding the total 
incident energy Ec.m., the total angular momentum J and the entrance channel 
mass asymmetry tjq. At the end of the subsection we will discuss the diabatic 
effect as well as the effect of the shell corrections on the value of Pqn for 
different entrance channels t]q. For the present study, several reactions leading 
to the 256No CN will be selected. 

Mesh density. Fig. 15 shows the dependence of some observables on AR (bot- 
tom) and Arj (top) for an entrance channel mass asymmetry 7]q = 0.0 and 
standard values of the remaining parameters. The total incident energy is 
Ec.m. = 30 MeV, while the total angular momentum is J = h. 

The results are quite stable regarding the values of AR (bottom) and Arj 
(top). The larger the mesh density, the larger is the transition probability 
rate between two neighbouring configurations (nodes of the mesh) at fixed 
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Fig. 13. (Color online) Intrinsic excitation energy (in MeV) of the nuclear shapes in 
Fig. 11. See text for further details. 



t = 7*l0"''''s t = 20*10"'''' s 




Fig. 14. (Color online) The same as Fig. 13, but the nuclear shapes are those in Fig. 
12. See text for further details. 



34 




0.15 0.25 -1 -0.5 0.5 1 

AR (fm) Tl 



Fig. 15. (Color online) Dependence of PcN-, Tqf and the mass distribution of the 
quasi-fission fragments on AR (bottom) and Ar] (top) for an entrance channel mass 
asymmetry rjQ = 0.0. See text for further details. 

values of the Gaussian variances and a^. It is reflected in a slight increase 
of the quasi-fission probability at large fragmentations 1] and in a reduction 
of the quasi-fission time (bottom left). The large mass yield for r] very close 
to one (top right) is associated with the strong decreasing of the asymptotic 
adiabatic PES at very large rj (see Figs. 8 and 10). The decay of the system 
for large R clearly determines the quasi-fission time which correlates with the 
AR values. This is because the difference between the level density of two 
neighbouring configurations is much bigger in the quasi-fission region than in 
the other regions due to the Coulomb repulsion energy. The main peak in 
the quasi-fission mass distribution corresponds to the entrance channel mass 
asymmetry r/o- The structures in the mass yield are related to the structures 
of the dynamical PES (see Fig. 10). These latter are initially caused by the 
diabatic effects, and by the remaining shell effects after the relaxation of the 
diabatic PES to the adiabatic one. The shoulders in the quasi-fission yield 
around t] ~ ±0.4 are related to the fragmentation ™Ni + ^^^W, which refiects 
the importance of the light fragment closed proton shell Z = 28. This is also 
consistent with the experiments for ^^Ca + ^°^Pb discussed by Itkis et al. in 
Ref. [15]. 

Critical parameters. Fig. 16 shows the dependence of the same observables 
discussed above on the parameters Kq (bottom) and Fq ^ (top) for an entrance 
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Fig. 16. (Color online) Dependence of the same observables of Fig. 15 on kq (bottom) 
and Fq ^ (top) for the entrance channel mass asymmetry riQ = 0.0. See text for 
further details. 

channel mass asymmetry rjo = 0.0 with Ecm. = 30 MeV and J = h. The 
remaining parameters have standard values. 

At the top in Fig. 16 (left) we can see that Pcn is more sensitive to the value 
of of Fq ^ than to the value of kq (bottom left). The opposite happens with the 
quasi-fission time Tgf. However, taking into account that here the Pcn scale 
is large, Pcn remains rather stable in the two dependences. The smaller Fg ^, 
the slower the diabatic PES relaxes. Thus more probability moves into the 
quasi-fission region, decreasing the value of Pcn, before the repulsive core of 
the diabatic PES at small R disappears. In this case the quasi-fission mass 
distribution (red curve at the top right pannel) notably reflects the structures 
of the initial diabatic PES (see Fig. 10). By decreasing Fg ^ by a factor ten, the 
quasi-fission time Tgf increases by a factor two. It is because a significant prob- 
ability still remains in the competition region after the complete relaxation of 
the diabatic PES to the adiabatic one. 

Although Pcn very weakly depends on Kg (bottom left), the quasi-fission time 
Tgf as well as the quasi-fission mass yield (bottom right) are strongly affected 
by the value of kq. The quasi-fission time Tgf clearly correlates with the kq 
value. The larger kq, the larger is the transition probability rate between the 
configurations, particularly in the quasi-fission region, due to the level density 
(phase space) effect discussed above in the dependence on AR values. For 
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a given period of time (with the relaxation time of the diabatic PES being 
determined by the fixed Fq ^ value) more probability moves into the quasi- 
fission region with increasing kq, thereby reducing the value of Pcn- The 
quasi-fission mass distribution (bottom right) reveals that its width decreases 
with increasing ko- 

Dependence on Ec.m., J tiq. Fig. 17 shows the dependence of the same 
observables discussed above on the total incident energy Ecm. (bottom, J — 
h and rjo — 0.0), on the total angular momentum J (middle, Ecm. = 30 MeV 
and 770 = 0.0) and on the entrance channel mass asymmetry 770 (top, J — H 
and Ec.m. — 30 MeV). The remaining parameters have standard values. 

At the bottom in Fig. 17 (left) we can observe that Pqn reaches a maximum 
value around Ecm. ~ 30 MeV, just at the Coulomb barrier of the diabatic cap- 
ture potential (see Fig. 9, bottom left), while the quasi-fission time remains 
constant. At energies above the peak of the Pcn excitation function, Pcn 
slightly decreases with increasing Ec.m.- After the relaxation of the diabatic 
PES to the adiabatic one in a time about 5 * 10~^^ s, a significant probability 
still remains in the competition region. The distribution of this probability 
between the fusion and quasi-fission regions depends on the competition be- 
tween the phase space of the two regions, which seems to be regulated by the 
incident energy Ec.m.- With decreasing Ec.m. towards the capture barrier, the 
phase space of the fusion region seems to play a more important role and, 
therefore, Pcn increases slightly. At low incident energies below the diabatic 
capture barrier, the fusing system spends a certain period of time in a clas- 
sically forbidden region. Consequently, the phase space of the quasi-fission 
region clearly dominates and Pcn decreases notably. 

At the bottom on the right in Fig. 17 it is shown that the mass yield of 
the quasi-fission fragments is correlated with the Pcn excitation function. At 
the lowest incident energy, the mass yield shows structures that reflect those 
caused by the diabatic effects in the entrance diabatic PES (see Fig. 10). 
With increasing Ec.m. these structures in the mass yield gradually disappear, 
although very smooth shoulders can still be observed at E^m. = 80 MeV. 
The quasi-fission mass distribution becomes smoother with increasing E^m. 
refiecting the asymptotic adiabatic PES (sec Fig. 10). The remainig smooth 
structures in the mass distribution at the highest incident energy are due to 
some remaining effect of the initial diabatic phase of the reaction along with 
the remaining ground-state shell effect of the warm quasi-fission fragments. 

In the middle in Fig. 17 (left) we see that Pcn globally decreases with in- 
creasing J and the quasi-fission time Tqf remains constant. The values of Pcn 
very weakly depend on J up to J = 40 because the rotational contribu- 
tion to the PES for compact shapes is small. Since the rigid-body moment of 
inertia, which the system reaches at the contact radius, is much larger than 
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Fig. 17. (Color online) Dependence of the same observables of Fig. 15 on Ec 
(bottom), J (middle) and r/o (top). See text for further details. 



the moment of inertia for the orbital motion at distances around the radius 
of the Coulomb barrier, the centrifugal effects increase the Coulomb barrier 
and it becomes narrower with increasing J. As a consequence, Pcn first in- 
creases slightly due to the increase of the Coulomb barrier and then (around 
J = 20 h) starts decreasing because the narrowing of the barrier clearly favors 
quasi-fission. For J > 40 h, the rotational contribution to the PES for com- 
pact shapes is so significant that these shapes are in a classically forbidden 
region, even after the complete relaxation of the initial diabatic PES. Thus, 
the dominance of the phase space of the quasi-fission region over the phase 
space of the fusion region is very strong. Hence, Pcn decreases strongly. The 
mass yield of the quasi-fission fragments (in the middle on the right in Fig. 
17) increases with increasing J. 

At the top in Fig. 17 (left) it is shown that Pcn is much smaller for near- 
symmetric reactions compared to very asymmetric ones. This is because the 
quasi-fission is stronger for more symmetric entrance channels due to the effect 
of the diabaticity that decreases the capture valley (see Fig. 9). The capture 
valley becomes deeper containing more compact shapes with increasing rjo. On 
the other hand, the quasi-fission time Tqf globally decreases with increasing 
?7o. Some structures are revealed, which are clearly correlated with those of 
the Pcn values (see Fig. 18). In contrast to the quasi-fission mass distribution 
for Tjo = 0.0 (black curve at the top on the right panel in Fig. 17), the mass 
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Fig. 18. (Color online) The same for Pqn as that in the top left panel of Fig. 17 
(black solid curve), but compared to values obtained without shell corrections (blue 
dashed curve) and without diabatic effects (red solid curve). See text for further 
details. 

yield of quasi-fission for tjq = 0.625 (very close to the system ^^Ca + ^"^Pb) 
is notably asymmetric with respect to 77 = 0.0. This asymmetry seems to 
increase with increasing 7]q (comparing the red curve for tjq = 0.2 to the green 
one for tjq = 0.625). It is important to remember that this happens in the 
rotating body-fixed reference frame. Further work is needed to clarify whether 
this anisotropy can be observed experimentally. The mass distribution for 
rjo = 0.625 also shows a local minimum instead of a peak at rjo, while Pcn 
reveals a local maximum (within 99% of accuracy of the calculation) in the 
small picture inserted in Fig. 18. Thus a reaction induced by closed shell nuclei 
seems to increase Pcn- 



Diabatic and shell corrections effect. In Fig. 18 (similar to the top left panel of 
Fig. 17) we can observe that the shell structure of the nuclei in the entrance 
channel plays an important role in establishing the value of Pcn- Here, Pcn 
as a function of r/o is shown for three types of calculation: (i) all shell effects 
(diabatic and shell corrections effect) are included (black solid curve), (ii) 
without shell corrections (blue dashed curve), and (iii) without diabatic effects 
(red solid curve). The effect of the diabaticity on Pcn can be seen comparing 
(i) and (iii), whereas the effect of the shell corrections comparing (i) and (ii). 
From this study we conclude: 
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• The effect of the diabaticity in suppressing Pcn can be very strong in near- 
symmetric colUsions. 

• The shell corrections can play a very important role in establishing the value 
oi Pcn- Depending on the entrance channel mass asymmetry, they can favor 
or inhibit the CN formation. 

The reaction with rjo = 0.0 (-"^^^Sb + ^^'^Sb), although more symmetric than 
that with rjo = 0.2, reveals a Pcn value about 0.96 due to the favorable 
effect of the fragment proton number that is very near the closed proton shell 
Z = 50. This value of Pcn is similar to those for very asymmetric entrance 
channels. However, the symmetric entrance channel shows a large diabatic 
shift of the capture barrier in Fig. 9 (bottom left panel) which inhibits the 
capture probability. Therefore, we can expect a very asymmetric entrance 
channel involving closed shell nuclei to be the best suited to form a heavy 
compound system at a given excitation energy. 



4 Concluding renicirks and outlook 

A theoretical formulation of the competition between fusion and quasi-fission 
in a heavy fusing system has been presented. Fusion and quasi-fission result 
from a diffusion process in an ensemble of nuclear shapes, each of which evolves 
towards the thermal equilibrium. The theory is unique and more realistic than 
the current fusion theories because it is based on microscopical grounds and 
also for the first time incorporates a wide range of important physical effects 
that impact on the formation of the heaviest compound nuclei. The dynamical 
collective PES emerges as a very important new feature of this approach, which 
partially reconciles conflicting aspects of the current models for CN formation. 

Realistic calculations for several reactions leading to ^^^No indicate that (1) 
the diabatic effects are very important in the capture stage of the reaction 
as well as in the subsequent evolution of the compact nuclear shapes in near- 
symmetric collisions, (2) these effects cause structures in the dynamical col- 
lective PES that are reflected in observables like the mass distribution of the 
quasi-flssion fragments and their intrinsic excitation energy. These structures 
are not related to the well-known static ground-state shell corrections, but to 
the diabatic sp motion through the shell structure of the nuclear shapes, (3) 
the quasi-flssion time is notably larger than the relaxation time of the initial 
diabatic collective PES to the adiabatic one, so the nuclear shapes reach the 
thermal equilibrium during the CN formation, (4) the compact fused shapes 
are much warmer than the quasi-flssion fragments, (5) the remaining shell cor- 
rections to the adiabatic PES can play a very important role in establishing 
the Pcn value, (6) Pcn very weakly depends on the total angular momentum 
J for small values of J (e.g., up to J ~ AOh for a central collision at the 
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capture barrier energy) and reaches its maximum value around the capture 
barrier energy, (7) the quasi-fission mass yield increases with increasing J, 
(8) the quasi-fission time very weakly depends on the incident energy and the 
total angular momentum, while it globally decreases with increasing entrance 
channel mass asymmetry. Here, some structures depending on the entrance 
channel mass asymmetry are manifested, which are correlated with those of 
the PcN values, (9) the angular anisotropy of the quasi-fission fragments may 
increase with increasing asymmetry of the entrance channel, and (10) very 
asymmetric reactions induced by closed shell nuclei seem to be the best suited 
to synthesize the heaviest compound systems. Last but not least, the calcula- 
tions are quite stable regarding their dependence on model parameters. 

The effect of other important collective degrees of freedom like the deformation 
of the colhding nuclei and their mutual orientations as well as the effect of the 
fission valley on the CN formation will also be investigated in the future. 
Works in this direction along with the calculation of ER cross sections are in 
progress. 
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A Derivation of the generalized master equation (13) 

The following coupled equations are obtained by acting with the projection 
operators C and Q = 1 — C on the Liouville equation (8) [|p(t)) = C\p{t)) + 
Q\p{t)) and replacing the partial derivative by a total derivative] 



CL{t)[C\p{t)) + Q\pm, 



(A.l) 



QL{t)[C\p{t)) + Q\pm. 



(A.2) 



Integrating (A.2) the solution Q\p{t)) reads 



Q\p{t)) = exp[-ig I dt'L{t')] ■ {Q\p{to)) 



to 

t tl 




(A.3) 
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Inserting (A. 3) into (A.l) and taking into account that CL{t)C = [51,52] 
because La^ult) = according to the definition (9), the diagonal part of the 
density matrix C\p{t)) obeys the following equation 



d 



C\p{t)) = -iCL{t)exp[-iQ J dt'L{t')]-Q\p{to)) 



dt 

to 

t-to 



-CL{t) j dTQ-x^[-iQ j dt'L{t')]QL{t - T)C\p{t - t)), 

(A.4) 



t-T 



where t = t — ti. The projection on {nn\ and summation over neH^,, taking 
into account expressions (10) and (12), leads to 

dP (t) 

^^ = Y: j dTK,{t,T)d,P,{t-T)+G,{t,to), (A.5) 
M 



where K^n{t,T) and G^{t,to) are given by expressions (14) and (15), respec- 
tively. Owing to the validity of the following relation 

J2 d^,K,t,{t, r) = ^ d,K,,{t, t)^-^^ d^K,^{t, r), (A.6) 



because 



E E ^»&.mm = 0, (A.7) 



according to definition (9), expression (A.5) turns into (13) after using the 
relation (A.6) in (A.5). The validity of (A.6) can be checked applying the 
multiplication rule for tetradics [51,52] regarding the product of the terms of 
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